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Abstract 
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■ Using the Calogero-Moser model and the Nakamura equations for 

i a multi-partite quantum system, we prove an inequality between the 

^1^1 mean bi-partite entanglement rate of change under the variation of 

a critical parameter and the level-curvature. This provides an upper 
bound for the rate of production or distraction of entanglement in- 
^ I duced dynamically. We then investigate the dependence of the upper 

C7^' bound on the degree of chaos of the system, which in turn, through 

the inequality, gives a measure of the stability of the entangled state. 
^ I Our analytical results are supported with extensive numerical calcu- 

V-j ' lations. 
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1. INTRODUCTION 



To store , process or teleport quantum information we need real, physi- 
cal quantum systems which live outside the "protected environment " of the 
quantum mechanics text books, and beyond the safe bounds of linearity 
These Quantum Information Devises (Q.I.D.) are always open quantum sys- 
tems and complex enough that the effects of non-linearity cannot be ignored. 
Given the fact that the prime and most relevant physical issue for QID is 
the support of quantum information, the main interest of most analytical, 
numerical and experimental work currently done is about the production 
and protection of entanglement, which is the manifestation of non-classical 
information in quantum systems. Since the physical devises are open sys- 
tems and may operate in a non-linear regime, it is a natural question to be 
asked whether the dynamics of production or distraction of entanglement is 
related to the dynamical properties of the corresponding classical systems. 
There are many works dealing with this problem, namely , the connection 
between entanglement and quantum chaos, El El El IH El El UDl 1121 US]- 
While it seems that there are some controversial results about the distrac- 
tive or constructive influence of classical dynamics on entanglement , one can 
classify and understand the conditions for one or the other behavior. The 
main problem with all these results is that they are model and initial state 
dependent. A better approach would use the Random Matrix ideas to ob- 
tain general results jTI]. Thus a general result in terms of inequalities and 
bounds can provide an insight for the relevant effects. Here we prove such 
an inequality which bounds the rate of change of mean bi-partite entangle- 
ment of a multi-partite system when a critical parameter is changed. The 
bound is proportional to the square root of level-curvature of the particular 
state of the system.Then using the level-curvature distribution function for 
different ensembles we find some optimality conditions which we believe are 
relevant for the design of QID. We made extensive numerical tests of the 
inequality using three different models. The inequality is satisfied strictly. 
Now , the established inequality can be applied to investigate the stability 
properties of the entangled state w.r.t. the degree of chaos in the system. 
The latter is quantified with the dependence of the upper bound on a critical 
parameter which controls the variation of the system from the regular to the 
chaotic regimes. We do that for two physically different systems. The first is 
described by a Hamiltonian which is the sum of a random matix with Pois- 
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son level spacing distribution and a random matrix with GOE distribution, 
[221 • The second system is a one dimensional spin chain with a random dis- 
tribution of defects [21]. First we confirm the transition from regularity to 
chaos by plotting the dependence of the parameter characterizing the degree 
of chaos vs the critical parameter of the dynamics. Then we compute the 
dependence of the bound on both the chaos parameter and the dynamics 
parameter. It is obvious from these plots , and the inequality , that the vari- 
ation of the mean bi-partite entanglement is bounded more strictly when the 
investigated systems are more regular. We believe that this is an indication 
of the usefulness of the inequality which has a model independent validity. 
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2. CHANGE OF ENTANGLEMENT AND AVOIDED CROSSINGS. 



We start from the autonomous non-integrable system with Hamiltonian 

H = Ho + tV (1) 

where Hq is integrable and tV a non-integrable pertubation with r the per- 
turbation parameter. Let the eigenvalue equation for each r be discrete and 
non-degenerate 

HiT)\niT) >= e„(r)|n(r) > (2) 

Then considering the parameter r as a 'time' the dynamical system (enir) , pn^r)) 
satisfies the equations [131 CHI UHl UZI 

=P„ (3) 



dr 

d\n{T) > 
dr 

d < n(T)\ 



dr 



2 E Vnm{r)V„^n{r){en{r) - eUr))-' (4) 
-- liTT'i^) > V.mnir){en{r) - em{r))~'^ (5) 

m(^n) 



dr 



- < m{r)\Vmn{r){en{r) - eUr))-' (6) 

where p„(r) = Vnn = < n{r)\V\n{T) >. These equations can be transformed 
into an integrable Hamiltonian system, the Calogero-Moser N-particle system 
in H-l dimensions. Now we suppose that we have a bi-partite system with 
subsystems A and B and the total system dynamics given by equation (1). 
Let the pure state 

Plsir) = |n(r)><n(r)| (7) 

Using the linear entropy El as a measure of entanglement we analyze its 
dependence on the parameter r. We have 

El = 1- tr{{plf ) (8) 

where 

Pa = trs^pls) (9) 



3 



We want to compute the derivative dEL/dr. We have 

§^tr[{pl{r)r] = 2tr[trB{p'XB)trB[^{plB)]] (10) 

Now 

< k{r)\^pUr)\l{r) > 

^< k{T)\{^\n{r) >) < n{T)\l{r) > + < k{r)\n{r) > < n{T)\)\l{T) > 

d d 
= < /c(r)|(— |n(r) >)5ni + 5kn{^ < n(r)|)|/(r) > 

and this , using Nakamura's equations above, is equal to 

7^,(r) = Snlil - dnk)-^ + 5nl{l - 4n)-3- (H) 

Thus 

i-plBir) = E7."IMT)></(r)| (12) 

kl 

and finally 

= -2j:rkitr[trB{\n{r) >< n{r)\)trB{\k{r) >< /(r)|)] (13) 

kl 

or 

= -2Y.rkMr) (14) 

with 

AI,{t) = ir[irs(|n(T) >< n{T)\)trB{\k{T) >< 1{t)\)] (15) 

We observe here that, while ^^/(r) is a bounded quantity, it can be proved 
that has first order poles at the avoided crossings, a fact we avoid assuming 
non-degeneracy of the spectrum. But if the system is classically chaotic 
we have many avoided crossings, and this leads to the result that for the 
corresponding states, at these critical values of the parameter r we expect 
fast change of the entaglement under a small change of the parameter. 
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3. CHANGE OF ENTANGLEMENT AND LEVEL-CURVATURE. 



Here we are going to show that the rate of entanglement change at the 
particular value of the parameter r is proportional to the square root of the 
level-curvature, under the assumption that for this value an avoided crossing 
of two levels gives so a small energy gap that only these two levels contribute 
to the sums above. We call these levels and 1. First we write equation (14) 
as 

From Nakamura's equations we have for the level-curvature 

- = i: (17) 

Under the assumption stated above we have 

dr €o — €i 

and similarly 

dEl _ Vo^Al, + V^oA\, 



dr 61 — 60 

For the curvatures we have 



2 ' 1" (19) 



|2 

6o — 6i 



Ko = 2 < (20) 



and 



|2 



^ 2^^ > (21) 



ei - eo 

that is under the assumption we have Kq — —Ki. We then obtain 



— ^ '. '. 

dT iHil' 
and eliminating the potential we get 

dEl 2V2Re[A%e''f'] 



Ko (22) 
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Ko\ (23) 



and similarly for the level 1. is a phase related to the potential matrix 
elements. Wc may note here that from the last two equations we have near 
the avoided crossings 

AEi ^ Re[Al,e^f] 

AEl Re[A%ei't'] ^ > 

a relation between the entanglement changes of the two levels. 



4. GENERALIZATION TO AN N-QUBIT SYSTEM. 



For an N-qubit system we use as a measure the mean bi-partite entangle- 
ment for the level n given by 



2 ^ 



(25) 



where pj = trALL^j{\n >< n]) Motivated by the above results one may 
conjecture that there is a similar relation between the rate of change of the 
average entanglement and some functional of level-curvature. It turns out 
that it is straightforward to obtain an inequality for the ground state in the 
form 



dr 



(26) 



Following the same procedure as that for the bi-partite entanglement rate we 
get the expression 



dr 



N 



(27) 

fe=i - eo 

where here we have —N < < +N We prove the inequality by taking the 
maximum of the r.h.s of (26). We have 



max- 



dr 



^ Re\Vok\ 

— 

k=\ - eo 



< 



2^-1 



k=i - eo 



2^-1 



^ £ ^fc - eo 



< 



2^-1 



Wc 



\ £ efc - eo i£l efc - eo 



2^-1 

E 



hx\Ka 



(28) 
(29) 
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where 



2^-1 

E 

k=l 



— ^0 



(30) 



Our numerical investigation, using Random Matrix Theory, indicated that 
this value of b is not optimal for the inequality. We may obtain a smaller 
bound under certain statistical assumptions, but still without involving par- 
ticular properties of the dynamics. The bound is , at least for our numerical 
results , amazingly good. The derivation is as follows. In (26) A^i is a sum 
of N terms for the N-qubit system 



JV 



^kl 



kl 



(31) 



Now we assume that these N terms are independent random varables uni- 
formly distributed in the interval [-a, a] with zero mean. From this it follows 
for the standard deviation cr{A^i) ~ Our assumption is to use for the 

derivation of the inequality the plausible bounds 



This lead to the expression 



h' = 



8a 



/3]V\ 



E - 

fe=i 



(32) 



(33) 



Now it is natural to suppose that a — l/N, but we get very good results 
taking the value a — 1/2^. We do not have an obvious answer for that. 



5. ENTANGLEMENT VARIATION AND THE DEGREE OF CHAOS OF 
THE SYSTEM. 



5.1 Dependence of the bound on the degree of chaos. 



Inequality (26) provides a general bound for the rate of change of mean 
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bi-partite entanglement. One may view this inequality as a criterion for the 
stability of entanglement under small changes of a critical parameter . The 
smaller the upper bound the more restricted the change of the entanglement. 
On the other hand a bigger upper bound cannot exclude small variation of 
the entanglement . For this we would need a lower bound. But still, it 
is interesting to investigate the dependence of the bound on the degree of 
chaos. Namely, it would be of some value to see whether integrability or 
chaos restricts more or less the possible change of entanglement or better 
whether there is an increased stability of entanglement which depends on 
the dynamical properties of the system. We use as a characterization of the 
dynamics the level statistic parameter j231 

/[P(s) - PwD]ds 

1 = i (34) 

J[Pp{s)-PwD]ds 



where sq ~ 0.472, which is the common point of Poisson Pp{s) = exp{—s) 
and Wigner-Dyson Pwd{s) = ^exp(— 7rs^/4) distributions, and P(s) the 
level density sistribution of the system. This parameter provides a measure 
of the degree of chaos, or better of relative distance of the level density from 
either the Poisson or the Wigner-Dyson distributions. We use two physically 
different models to analyze this question. The first is given by a Hamiltonian 
which is the sum of two random matrices, one with Poisson level statistics 
and the other with GOE, |23]. This superposition is a rotation which depends 
on an angle 6 which can be varied to give a continuous set of Hamiltonians, 
from a pure integrable to a pure chaotic one, in the RMT framework. We first 
study the dependence of 7 on ^ to establish the variation of the behavior of 
the system from the regular to the chaotic regimes. Then we make the plots 
of the dependence of the upper bound parameter b on both 7 and 9. These 
plots give a very prominent dependence of the bound on the dynamics of the 
system. We comment on this in more detail bellow. The second system is 
a one- dimensional chain of pair-wise coupled spins perturbed by a random 
distribution of defect sites [21]. These defects are introduced as differences 
in the two-level energies on these sites. Here the critical parameter for the 
transistion from the regular to the chaotic regime is the two-point correlation 
parameter d of the normally distributed random defects. We plot again the 
dependence of 7 on d to establish the behavior of the system. Then we plot 
the dependence of b, the parameter of the bound, on both 7 and d. We 
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observe the same behavior of the bound as with the first system. Again we 
comment in more detail bellow. 



5.2 Entanglement variation for different ensembles. 



We may apply the derived inequality to investigate the sensitivity of en- 
tanglement under small changes of the parameter r due to appriximations 
or environmental changes. The maximum change of the entanglement would 
be 

AQ = I^IAr = bJKoAr (35) 
or ^ 

Using this we may observe the differences in the influence between various 
symmetry ensembles. To do this we calculate the average of the entanglement 
change using the distribution functions of the level-curvatures, that is we 

compute [iHi cni Em EH 



We have 



< AQ > = b< ^jKo > At 

< ^/\K\>GOE^ 

< ^[\K\ >GUE ~ \/272 



(36) 

(37) 
(38) 



< y/\K\>GSE^ O.Q^^ (39) 

where 

7. = i^A (40) 

where A is an appropriate constant related to the ensemble density pT] . 
Then assuming a common value for b, this gives an estimate of the relative 
merit of the ensembles 

« 0.84 , « 0.7 (41) 

^WGUE ^WGSE 

Of course , both factors in the r.h.s. of the inequality depend on the dynam- 
ics, as it was stressed above, but still this result could be used as an indication 
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of the relative merit for different ensembles. In the designing strategy for de- 
vices which employ the property of entanglement the most serious problem 
is the fragiliry of resource. We believe that qualitative properties as the one 
in (41) could be a guide in this area. We currently investigate this point. 



6. MODELS AND NUMERICAL RESULTS. 
6.1 Testing the inequahty. 

We have checked the inequality (26) numerically choosing three different 
categories of models. All are of the form of (1) . 

Model A 

We take here 3 qubits with pairwise interactions and perturbed globally 
by a potential V. For Hq we choose 

3 

Ho = Yl (^j^^j + (42) 

j=l i<j 

where aj and A are arbitrary constants. Our choice gives b=9. 849476 and 
b'=0.410395. V is Hermitian with random matrix elements. 

Model B 

Hq is an arbitrary Hermitian matrix and V as in Model A. For 2 qubits we 
have b=11.66033 , b'=l. 190078 and for 3 qubits b=20.57092 , b'=0.8571215 
and for 6 qubits b=38.85935 , b'=0.1431131. 

Model C 

2 qubits with Hq as in Model B. Two choices for V, from GUE and GOE. 

For all models we computed the two sides of the inequality (26), for the 
given choice of Hq and for 3000 choices of V. For each pair of values we put 
a point in the graph of |^-| versus ^ \ Kq \ . In the same graph we draw the 

line y — b'J\Ko\. For all models we observe that all points, but for few, are 
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bellow this line. These outliers are due to the fact that we have estimated 
b' statistically. To make the results more obvious, for each case we draw a 
second graph for the difference of the l.h.s minus r.h.s of inequality (25). i.e. 

S = \^^\numericai — b^J\Ko\ . We Call 6 'saturatiou index' since it indicates 
how far we are from the saturation of the inquality. 




Figure 1: Model A: 3 qubits. Numerical test of (26). Left : distribution of 
3000 random points. Right: dependence of the saturation index on curvature 
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Figure 2: Model B: 2 qubits. Numerical test of (26). Left : distribution of 
3000 random points. Right: dependence of the saturation index on curvature 




Figure 3: Model B: 3 qubits. Numerical test of (26). Left : distribution of 
3000 random points. Right: dependence of the saturation index on curvature 
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Figure 4: Model B: 6 qubits. Numerical test of (26). Left : distribution of 
3000 random points. Right: dependence of the saturation index on curvature 




Figure 5: Model C: 2 qubits. Numerical test of (26). Left : distribution of 
3000 random points. Right: dependence of the saturation index on curvature. 
GUE perturbation 
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Figure 6: Model C: 2 qubits. Numerical test of (26). Left : distribution of 
3000 random points. Right: dependence of the saturation index on curvature. 
GOE perturbation. 

These results seem to confirm that the inequality is satisfied strictly. The 
plots of |^-| vs \Ko\ for all cases make it obvious that the inequality is cer- 
tainly not saturated, but the plots of 6 (the saturation index) vs \Ko\ gives 
us an indication how this lack of saturation depends on dynamics. It is in- 
teresting that the higher the value of \Ko\ the stricter the bounding of the 
entanglement change is. Of course for any specific model, both factors of the 



r.h.s of the inequality , namely b and ^J\Ko\ must be taken into account, and 
for comparison of the models we must , somehow, normalize this dependence. 
We report elsewhere on this point. 

6.2 Dependence of the upper bound on the degree of chaos. 

We used two models to investigate how the parameter b of the upper 
bound of the proved inequality depends on the dynamical properties of the 
system with bi-partite entaglement. These models have different physical 
interpretation, one being a general RMT picture of chaotic systems, and the 
other a simplified spin-chain with defects. As it turned out, for both sys- 
tems, the upper bound parametr b strongly depends on the degree of chaos 
in the dynamics of the system. This degree is quantified by the widely used 
parameter 7. Certainly more measures should be used, but we feel that , as a 




14 



first step, the present results are very characteristic. We present the models, 
some numerical details and the results with some comments. 

Model D 

The Hamiltonian of the model is the weighted sum of two random matri- 
ces, one with Poisson level statistics and the other with GOE [221 ■ 

H{e) = cos{e)Hp + sm{e)Hw (43) 

We have chosen the dimension of the matrices to be 2^ x 2'' to correspond to 
a system of 7 qubits. For ^ = we have regular behavior and for 6' = 7r/2 we 
have chaos. Thus we may vary 6 from to 7r/2 and check through the level 
spacing distribution parameter 7 the variation from the regular to the chaotic 
regimes. For each value of 6 we used a 100 different random realizations of 
the Hamiltonian and calculated the averages < b > and < 7 >. We took a 
special care of the data to deal with the outliers, applying a small variation 
of the methods in |24j. To compute 7 we had to fit the data appropriately. 
We have chosen to do that using the Weibull distribution. 

ya,c{x) = acx^^^exp{—ax^) , x G [0, +00] (44) 

The Weibull distribution contains as a special case the Brody distribution, 
which is usually applied to interpolate between the Poisson and the Wigner 
distributions. We use this model to see the dependence of the bound b on 
the degree of chaos. The model does not offer us the posibility to calculate 
the mean bi-partite entanglement, but it gives us the general behavior of b 
w.r.t. the underlying dynamics. In Figures 7, 8 and 9 we plot respectively 
the dependence of 7 and d on 6', and of b on 7. The first gives the picture of 
the change of the behavior of the system from regularity, high 7 , to chaotic- 
ity , small 7. The other plots give the dependence of the bound on degree of 
regularity or chaoticity of the system. It is obvious that b takes its smallest 
values in the region of regularity. This means that the inequality is bounded 
more strictly, that is, that the entanglement change is more restricted when 
the system is more regular. 



Model E 
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Figure 7: Model D: 7 qubits. Dependence of the level statistics parameter 
7 on the dynamical parameter 9. The diagram shows the variation from 
regularity to chaos as a function of the critical parameter. 



28 n 




18- 



-0,2 0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 

9 

Figure 8: Model D: 7 qubits. Dependence of the bound parameter b on the 
dynamical parameter 9. The diagram shows the variation of the bound of 
the inequality as a function of the critical parameter. 
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Figure 9: Model D: 7 qubits. Dependence of the bound parameter b on the 
level statistics parameter 7 . The diagram shows the variation of the bound 
of the inequality as a function of the measure of chaos. 

We use this model to test the dependence on the dynamics of both the 
bound parameter b and the mean bi-partite entanglement. This gives us 
the picture of how strict the inequality becomes and in the same time how 
much the entanglement is influenced by the underlying dynamics of the sys- 
tem. This model is well studied in the literature, |21j. It describes a one- 
dimensional spin chain with nearest neighboor interaction with strength J. 
These qubits are coupled to an external homogeneous magnetic field. The 
system is perturbed by a distribution of random defects in the form of random 
interaction energies hj with the external field, with normal distribution 

< hj > = , < hihj > = (f5ij (45) 

N J N-1 

H = Y,{h + h^)cj,, + (46) 
i=i ^ i=i 

It is known [24J, that the system is chaotic for the the values of d in the 
interval 0.1 < d < 0.5 for J=l. We have analyzed this system for the 
interval < d < 2.5 for a chain of 9 qbits. For each value of d we used 100 
random matrices to compute the averages < Q >, < 6 > and < 7 >, again 
using the WeibuU distribution to fit level distribution density. The results 
are presented in the Figures 10,11 and 12. In the two parts of Figure 10 
we plot 7 and Q vs the critical parameter d. In Figure 11 we plot b vs d 
and 7, and in Figure 12 we plot Q vs 7 and b. Figure 10 exhibits how the 
dynamical parameter d influences both the appearance of chaotic behavior 
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and the establishing of entanglement. Figure 11 shows the dependence of the 
bound parameter b on the critical dynamical parameter and on the measure 
of chaos in the system, and Figure 12 gives the dependence of mean bi- 
partite entanglement on the measure of chaos and on the bound of the rate 
of change. From Figure 10 we see that at about the value of d=0.25, 7 is 
minimal and Q is maximal. This means that at this value we are closest to the 
chaotic behavior while in the mean time the mean bi-partite entanglement 
is maximal, something that conforms to similar results in the literature. 
From Figure 11 it follows that at this value of d the bound parameter b 
is maximal, something which means, that while the entanglement becomes 
maximal, the inequality bound parameter becomes less strict. Or, in other 
words, that the entanglement is more fragile. Finally Figure 12 reveals the 
direct relation between entanglement and the bound parameter b. It shows 
that more entanglement is associated to less strict bound , or to a more fragile 
regime. Of course, in order to draw conlusive results, we need to use bigger 
matrices and certainly more realistic systems. But , even at this stage, our 
results point to the usefuUness of our inequality for giving a general picture 
of the situations where we have the interplay of classical chaotic dynamics 
and the quantum regime of entanglement. 




18 



100 n 






















100-, 


90- 






















90- 


80- 

70- 






















80- 
70- 


60- 






















eo- 


50 






















K 50- 



0,0 0,5 1,0 1,5 2,0 2,5 0,5 0,6 0,7 0,8 0,9 1,0 

d y 

Figure 11: Model E: 9 qubits. Dependence of the bound parameter b on the 
dynamical parameter d and on the measure of chaos 
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Figure 12: Model E: 9 qubits. Dependence of the mean bi-partite entangle- 
ment on the measure of chaos and on the bound parameter b 



7. DISCUSSION. 



We have established a level-curvature dependent bound for the change 
of mean bi-partite entanglement under small changes of perturbation of the 
system. Extensive numerical tests were performed indicating no violation 
of this bound for a selected class of models. This bound depends on two 
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factors, the parameter b and the curvature of the levels. Both factors reflect 
general properties of the undelying dynamics and combine to bound the rate 
of change of the mean bi-partite entanglement under small variations of a 
critical dynamical parameter. Using five models we made extensive numerical 
investigation to check various aspects of the proposed inequality. First we 
tested the inequality with three different models and we found no violation. 
The plots of |^-| versus l^ol indicate that though the inequality is satisfied, 
it is far from being tight. To test the dependence of this lack of tightness 
, we plotted a saturation index vs \Kq\. From these plots it is evident that 
the saturation is better with increasing mean curvature. The second type of 
tests concern the dependence of the bound parameter b on the dynamical 
property of regularity or chaos of the system. For two different models we 
first checked the dependence of the level statistic parameter 7 on the critical 
parameter of the models. We use 7 as a measure of the distance of the model 
from being integrable or chaotic. Then we plotted the parameter b of the 
bound vs both 7 and the critical dynamical parameter. It is evident from 
these plots that there is a strong dependence of the bound on the degree of 
chaos. Since b is computed for the unperturbed system, the size of the bound 
reflects somehow the fragility of entanglement under small perturbation of 
the values of the dynamical parameters. Finally the plots of the dependence 
of entanglement vs the measure of chaos and the bound parameter, give a 
first picture of how the proposed inequality could be used to approach the 
fragility of entanglement. The models we analyzed are general enough, and 
to our opinion, cover the properties of many specific models used in the 
literature. Our results, we believe, would be very useful for the study of 
the influence of inaccuracies of perturbation parameters or the influence of 
stochastic environments. Furthermore, using the curvature, one could design 
optimal strategies for entanglement and protection from decoherence. We 
currently investigate the use of the bound for models of concrete realization 
of Q.I.D.'s. 
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